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Abstract 

We study the cosmological evolution for a universe in the presence of a continuous tower of massive scalar 
fields which can drive the current phase of accelerated expansion of the universe and, in addition, can 
contribute as a dark matter component. The tower consists of a continuous set of massive scalar fields with 
a gaussian mass distribution. We show that, in a certain region of the parameter space, the heavy modes 
of the tower (those with masses much larger than the Hubble expansion rate) dominate at early times and 
make the tower behave like the usual single scalar field whose coherent oscillations around the minimum 
of the potential give a matter-like contribution. On the other hand, at late times, the light modes (those 
with masses much smaller than the Hubble expansion rate) overcome the energy density of the tower and 
they behave like a perfect fiuid with equation of state ranging from to -1, depending on the spectral index 
of the initial spectrum. This is a distinctive feature of the tower with respect to the case of quintessence 
fields, since a massive scalar field can only give acceleration with equation of state close to -1. Such unique 
property is the result of a synergy effect between the different mass modes. Interestingly, we find that, for 
some choices of the spectral index, the tower tracks the matter component at high redshifts (or it can even 
play the role of the dark matter) and eventually become the dominant component of the universe and give 
rise to an accelerated expansion. 



1. Introduction 

More than a decade after the first evidences for the accelerated expansion of the universe [l[ , the true 
underlying mechanism responsible for it still remains unclear. After all this time, a simple cosmological 
constant still seems to succeed in explaining most of the cosmological observations by using one single 
additional parameter 0]. This has lead to establish the so-called ACDM model as the standard model of 
cosmology, in which we include the cosmological constant A accounting for about 70% of the energy content 
of the universe and the Cold Dark Matter component comprising 25% of the universe content, that amounts 
to 85% of the matter component. Despite its remarkable phenomenological success, alternatives to the 
standard ACDM model have been extensively explored because it is not completely satisfactory from a pure 
theoretical point of view. In particular, we need to assume that 95% of the cosmological budget is due to 
unknown fields or exotic forms of energy. Concerning the dark energy sector, alternatives to explain the 
accelerated expansion have been considered in part due to the tiny required value for A, which turns out to 
be many orders of magnitude smaller than its expected natural scale (see Q for an extensive review on the 
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cosmological constant problem). This naturalness problem is also linked to the coincidence problem referred 
to the fact that dark matter and dark energy give comparable contributions to the total energy density of 
the universe precisely today, while they have evolved very differently throughout the universe expansion 
history. 

A simple step forward is to assume that the accelerated expansion is not caused by a cosmological 
constant, but it is driven by some dynamical field Q whose evolution could solve the naturalness and/or 
coincidence problems of the cosmological constant. The most widely studied alternative models for dark 
energy are by far those based on scalar fields with a certain potential (quintessence) or non-standard kinetic 
terms (K-essence) very much like in the inflationary models for the early universe. However, also 
higher spin fields have been considered as dark energy candidates, like vector fields @ or general p-forms 
0. The difference between inflationary models and dark energy models are the involved scales and that, 
whereas during inflation the accelerated expansion must end, for dark energy models is desirable to have the 
accelerated solutions as attractors so that the natural evolution of the universe is towards an accelerated 
phase. Besides its simplicity, scalar fields are a common feature of high energy physics that extend the 
standard model of particles, so it is natural to think that they can be present in the universe and play a 
relevant role in the cosmological evolution. An alternative explanation for the cosmic acceleration could 
also be that the gravitational sector is modified on large scales, i.e., the validity of General Relativity could 
breakdown in the infrared regime. In this context we could mention the so-called f{R) theories Q or those 
based on extra dimensions like the DGP model Q. It is however a typical feature of modified theories of 
gravity that the new effects can be captured by a certain scalar degree of freedom like df{R)/dR in the 
f{R) theories or the bending mode in braneworld models. This is so because any modification of General 
Relativity necessary leads to the appearance of new degrees of freedom and this generally results in a scalar 
field. 

Although the kind of models mentioned in the previous paragraph can give equally consistent explana- 
tions for the cosmic acceleration as a cosmological constant does, naturalness problems usually still remain. 
For instance, in the simplest case of quintessence with a canonical massive scalar field, we need its mass to be 
m < Hq — 10"'^^ eV, again a tiny value when compared to what one would expect from high energy physics 
Q. Moreover, such a value is unstable when including quantum corrections. This constraint on the effective 
mass is necessary for the potential energy to dominate over the kinetic energy of the field today and thus 
its equation of state is near —1. The situation does not improve much by resorting to more sophisticated 
potentials or non-standard kinetic terms. Most of the times, the naturalness problems reappear hidden 
behind some parameters of the potential. 

On the other hand, dark matter models resorting to scalar fields also exist. For instance, the lightest 
massive mode of a Kaluza-Klein (KK) tower corresponding to the compactification of one extra dimension 
is a viable candidate for dark matter. Stability of such a particle is guaranteed by KK-parity corresponding 
to momentum conservation associated to spatial translation in the fifth dimension. Also resorting to extra- 
dimensions we can mention the vibrations of our brane in the extra dimension (branons) that can play the 



role of dark matter [13| . More recently, an ensemble of scalar particles which can interact among them has 



been considered as a dark matter candidate . The interesting feature of this approach is that the scalar 
particles do not need to be stable, but a certain balance between the width of the decays and the abundances 
allows to have a dynamical scenario in which the different effects conspire to give an overall dark matter 
contribution. In [l2|, it was also considered a scenario with multiple Kaluza-Klein dark matter candidates. 

In this work we present a novel approach to the dark energy problem based on scalar fields that, in 
addition, allows to unify it with the dark matter component. Unlike previous studies, we do not rely on one 
single scalar field, nor a finite discrete collection of (perhaps interacting) scalar fields. Here we shall study 
a cosmological scenario with a continuous tower of massive scalar fields which only interact with each other 
by means of gravity. Towers of scalar fields arise in a natural manner in nonlocal theories 3, 15, 1^ or 
theories with extra dimensions as we have discussed above. If one considers one compact extra dimension, 
a massless scalar field leads to a spectrum in the 4 large non-compact dimensions consisting of a massless 
scalar field plus a discrete tower of massive scalar fields, being the corresponding masses determined by the 
compactification scale of the extra dimension. Moreover, the mass gap between two consecutive modes of the 
tower is inversely proportional to such a scale. Thus, as the size of the extra dimension gets larger, the mass 
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gap between the different components of the tower gets smaller and. in the limit when the extra dimension 
is infinite, the tower becomes continuous. This actually happens in models with extra large dimensions like 
the Randall-Sundrum or DGP models. Another context in which a continuous distribution of massive fields 
appear is as an effective description of the unparticle fields [l^l • 

Here however we shall adopt a more phenomenological approach and shall simply assume the presence of 
a continuous tower of massive scalar fields in our universe, without paying special attention to the underlying 
mechanism generating it, that could be either unparticle physics or braneworld scenarios, or even more exotic 
set-ups. Our aim is to show that a continuous tower of simple massive scalar fields can lead to accelerated 
expansion and, in addition, it can simultaneously give a dark matter contribution and, thus, allowing for 
the unification of both dark components of the universe. 

The paper is organised as follows. In the next Section we give our starting action for the continuous 
tower that we shall use for our analysis and we will derive the relevant equations. In Section 3 we study 
the cosmological evolution of the tower during the standard radiation and matter dominated epochs and 
show how the light modes of the tower eventually dominate the universe and can produce an accelerated 
expansion. In Section 4 we study the case in which the tower dominates the universe and show how we can 
indeed have a transition from matter domination to accelerated expansion in a tower dominated universe. 
Finally, in Section 5 we discuss the results of the paper. 



2. A continuous tovifer 

As we have mentioned in the Introduction, we are not concerned about the true underlying mechanism 
generating the tower, but we shall adopt a phenomenological approach. Thus, we simply assume that our 
tower is described by the action 

£^ - i y dmV{m) {df,^rad^(j)m - rn^iQ (1) 

where V{m) is the spectral mass distribution of the tower, that we assume to be gaussian with zero mean 

V{m)^-^e-f^. (2) 

It is interesting to note that the only free parameter present in the action is the dispersion of the distribution 
a. Notice that, although m runs over the real numbers, no tachyonic modes are actually present because 
is always positive. In fact, since the mass distribution peaks at m = 0, the corresponding equations will 
remain invariant under the change m — >■ — m, i.e., we have a mass reversal symmetry and we could consider 
just m > without affecting the physics of the system. One could of course consider different spectral mass 
distributions like a gaussian distribution for m^, a x^-distribution or even more sophisticated distributions. 
However, the key point for our purposes here is that the distribution peaks at the origirf^ and has a given 
width, being the particular details of the precise distribution under consideration not really relevant. The 
reason for that is that we want the lightest modes of the tower to play a relevant role in the cosmological 
evolution so that they can eventually drive an accelerated expansion, as we show below. Actually, having 
a distribution such that the origin is within the 1-cr region would suffice for the very light modes not to be 
irrelevant. 

On the other hand, we are also considering the simplest case in which the different modes of the tower 
are decoupled from each other. In the most general set-up, the different mass modes could be coupled. 



^This also prevents the introduction of additional dimensionfuU parameters in the theory to specify where the distribution 



peaks. In our case, one could use P(r?i) oc exp 



(m-mp) 
2^ 



that depends on the mass mo giving the position of the maximum 



of the distribution. Even though we do not consider this possibility here (in order not to introduce unnecessary complications 
in our analysis), as long as mo > a our findings will still be valid. 
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although one can redefine the fields to use normal coordinates so that we end up with the above diagonal 
action. More explicitly, we could start from the action 

S ^ — f d^x\/—q f dm f dn 



Then, we can redefine the fields as follows: 

1pm = I A{m,p)(j)pdp (4) 

with A{m,p) satisfying 

dpdqA{m,p)A{n, q)K{p, q) = V{m)S{m — n), 

dpdqA{m,p)A{n,q)V{p,q) = m^V{m)5{rn — n) (5) 



with 5{m — n) the delta function. Then, it is straightforward to see that this redefinition diagonalizes our 
action so that we end up with 

5 = i j d^x^/^ j dmVim) [d ^(jymd^ (pm " ■rn^4>m4'm] (6) 

which is our initial action. Notice that this is, in principle, doable only for the quadratic part of the action. 
However, if we had interaction terms, the different modes will be unavoidable coupled. On the other hand, 
one could always redefine the fields to use canonically normalised fields (f> — !• 0/ ^JV{m) so that we get rid 
of the distribution V from the action. Then, the effect of the mass distribution can be transferred to the 
initial conditions. However, we prefer to keep V^m) in the action so that we can differentiate its effect form 
that of the initial conditions. 

In the present case, each mode satisfies a decoupled Klein-Gordon equation 

(□ + m2)0,„ = O. (7) 

In fiat spacetime, each component of the tower evolves like a free field and they are not affected by each 
other. The same happens in a curved background when they are treated as test fields, i.e., we ignore their 
backreaction on the spacetime. However, when their gravitational contribution is relevant, they become 
coupled by means of gravity through Einstein equations. 

The energy-momentum tensor of the tower is simply given by the superposition of the energy-momentum 
tensor of its components weighted by the mass distribution 

T^t = j dmV{m) [a^0™5,0„, - g^, {d,,dprad^(l,m - m'4>D ] ■ (8) 
Thus, the energy density and pressure are: 

dmV{m) {d^cPrndf^r^ + m^(t>l) (9) 
= J dmV{m) (9^<^™a^0"' - m^™) (10) 

so that the equation of state reads: 

_p^ _ Jdmrjm) {d^cPmd^f" - m^cPl) 

JdmV{m)id^cf>rnd^(j>"^+mm)- ^ ' 

Now that we have all the equations that we need, we proceed to study the cosmological evolution of the 
tower in the next Section. 



^Assuming that both operators K{m, n) and V{m, n) are symmetric it is always possible to diagonalize the action. 
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3. Cosmological evolution 

In a Friedmann-Lemaitre-Robertson- Walker (FLRW) universe described by the line element ds^ = dt^ — 
a{t)'^d3F, with a{t) the scale factor, the homogeneous mode of each component of the tower evolves according 
to the equation 

4>„, + 3H<j>,n + m^(f>m = (12) 

where H(t) — d/a is the Hubble expansion rate. Here we sec that the different modes only couple to each 
other when the tower contributes in a non-negligible way to the Hubble expansion rate. Since, in this Section, 
we shall assume that the universe expansion is driven by some background barotropic fluid (radiation or 
matter), the different modes will evolve independently of each other. 
The energy density in this case reads 

p^^^Jdm Vim) (,/)2„ + m^q^l^ . (13) 

If we now consider a power law expansion with H = p/t (as it happens during the radiation and matter eras 
with p = 1/2 and p = 2/3 respectively), we find the usual result 

(l)m{t) = t-" [A,M\mt\) + BM\mt\)] (14) 

with 1/ = (3p — l)/2, Ji, and the Bcssel functions of first and second kind, and Am and Bm integration 
constants that will, in general, depend on the mass of the corresponding mode. In particular, for the 
radiation era (p = 1/2) we have 

-^fn W = t-'^' [A^„J,/i{\mt\) + B^Y,^,i\mt\)] , (15) 

whereas for the matter era (p = 2/3) we find 

cb^^it) = {A'^l cos(|mt|) + Bll sin(|mt|)] . (16) 

The different masses of the modes will make them evolve in different ways so that they will contribute 
differently to the tower energy density. Depending on the ratio between the mass and the Hubble expansion 
rate, we have two regimes. For m ^ H (mt 1), (jjm is approximately constant, whereas for H <^ m 
{mt ^ 1) the field will oscillate with decaying amplitude. Thus, since the value of H decreases as the 
universe expands, the natural behaviour of a given component of the tower is to remain frozen (or, more 
precisely, in slow roll phase) at early times until H ^ m, when the field starts oscillating around the minimum 
of the quadratic potential. 

In the regime in which the field remains frozen, the energy density of the corresponding mode is constant, 
while in the oscillating regime the field behaves like a matter fiuid with an equations of state that averages 
to zero over one oscillation, i.e., w = 0. This is the well-known result in the usual quintessence models. 
However, having a continuous tower of scalar fields, this is not the end of the story, since the total energy 
density (or the effective equation of state) comprises the contribution from fields with any mass. In fact, 
this is a crucial difference because the problem for one single scalar field to play the role of dark energy is 
that its mass should be unnaturally small (~ 10"'^^ eV). For the case of the continuous tower, we have all 
the possible masses and, in particular, those tiny masses are always there. This is true for our distribution 
choice, because for distributions that peak at some finite mass mo and such that the origin is far away from 
the maximum, the light modes will be very suppressed and their role will be irrelevant. This is the reason 
why we need the peak not to be many sigmas away from the origin. 

According to the discussion of the previous paragraph, we can decompose the tower in two sets of modes, 
those with masses smaller than the expansion rate (that will evolve as a cosmological constant) and those 
with masses larger than the expansion rate (which will contribute as a matter fiuid). However, the overall 
behaviour of the tower will not be as simple as this picture might suggest, since we have additional effects. 
To see these effects more clearly, we will split the energy density of the tower as the contribution of the 
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light modes plus the contribution of the heavy modes, where by light and heavy we mean with respect to the 
expansion rate. In other words, we can write p = pi + ph with 



Pi ^ dmT{m)pra(t) and Ph ^ I dmV {m) pmit) . (17) 

J\rn\<H J\m\>H 

First of ah, only those modes with masses m < a will contribute in a non-negligible way to the tower 
because of the exponential suppression due to V{m). This means that ph is determined by the modes with 
H < m < a. In particular, if _ff > a, ph will be exponentially suppressed. However, since H decreases 
throughout the universe evolution, ph will acquire contributions that are not suppressed by V{m) after the 
Hubble expansion rate crosses the value of the tower dispersion a. Obviously, the cosmological evolution of 
the full tower depends on how the different mass modes contribute to the full energy density of the tower 
at each time and this will depend on the initial distribution. In order to set the initial conditions for the 
tower we should know the underlying mechanism leading to the generation of the tower so that we can 
compute the actual primordial spectrum of the tower. However, since we are not concerned here with the 
precise mechanism generating the tower, we shall assume a general power law primordial spectrum so that 
the initial conditions will be assumed to be given by 

a b 

(jjmiUni) = A — (18) 

m 

with A some amplitude and h the spectral index of the primordial tower distribution. Moreover, the field 
will be assumed to be initially at rest so that 4'm{tini) ~ 0. On the other hand, since we want both the light 
and the heavy modes to contribute in a non-negligible way at the initial time, we shall further assume that 
(7 > Hini t~^j^. Finally, we should note that, in order to avoid divergences in the energy density of the 
tower we must impose b < 3/2. To understand why this is so, we compute the initial energy density 

Pini = - j dm m^(f>l^{t„,,)V{m) -^-j= J dy-^^ (19) 

where y = The integral will not have convergence problems for y — > oo thanks to the exponential 
suppression due to V{m). However, we can find divergences near the origin because the integrand behaves 
like fj^g^f^ diverges for 2(6 — 1) > 1. Precisely for this reason, we shall restrict to the case with 

6 < 3/2. In FiglTJ we show the initial contribution of the different modes to the tower energy density for 
different values of b. 

It is clear that pt will always give a contribution decaying like a matter fluid, i.e., ph because, for 

our initial conditions which assume that a 3> Hini, the dominant contribution to ph comes from the modes 
with m ^ a. In fact, if we approximate the gaussian by a step function we have: 

A^a^ r dy A^a^ ^^^^ 



V2^ Jo y'^"-'^ V2^(3-26)- 

Then, the fact that H varies (and more modes contribute to ph) will not affect its cosmological evolution. 
Moreover, since we want the heavy modes to dominate over the light ones at the initial time (we shall make 
this statement more precisely below), we can write 

Ph ~ p^n^ ^ /— -— • (21) 



a 



/2^(3 - 2b) 



In this expression we see again why we must restrict to cases with b < 3/2. Otherwise, we would find a 
divergence. 

For the light modes however the variation of the amount of modes contributing to pi will impact its 
cosmological evolution because the corresponding integral will also be dominated by the upper limit which 
now is given by H{t), a time-dependent quantity. Physically, it is easy to understand the mechanism 
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Figure 1: In the left panel we show the initial energy density spectral distribution, i.e., 'P{m)pm{ti,ii) for b = 0.5, 1 and 1.2 
(long-dashed, solid and short-dashed respectively). There are two qualitatively distinctive initial distributions depending on 
b. If b < 1, the distribution peaks at some non-vanishing value, whereas for b > 1, the distribution peaks at the origin. In 
the right panel we show the evolution of the tower energy density in a matter dominated universe. We have chosen 6=1 
and a = lOHi„i so that the heavy modes that behave like a matter fluid dominate initially. We also show the evolution of 
individual components of the tower (dashed lines) to illustrate why the modes that are initially frozen and start decaying like 
a matter fluid at some time lead to a slower decay than the one corresponding to a single massive field. As we can see, the 
mechanism for this to happen is that we have a continuous tower with all the possible masses so that there will always be a 
lighter field that will decay later, making the overall tower energy density to have an equation of state different from 0. 



governing the evolution of the hght modes contribution. We show it expheitly in the right panel of FiglU 
where we plot the evolution of the entire tower in a matter dominated universe and also show the evolution 
of single modes. We see how the modes that are initially frozen when their masses are larger than H, give 
a matter-like contribution when H becomes smaller than the corresponding mass. However, since we have 
a continuous of mass modes, we have modes that start decaying at any given time. In other words, if a 
mode of mass m starts decaying at time tm, then there is a mode with mass m + Am that starts decaying 
at + A.tm so that the final contribution to the energy density is the envelope of all the individual energy 
densities, as shown in the figure. We can also estimate this analytically. To compute the evolution of the 
light modes that can give a dark energy contribution, we write 

Pi^t; dm mm^{U^i)V{m) ~ - = p„,, - (22) 

2 J|m|<H V27r(3 - 26) V cr y \a J 

where we have used that the light modes remain approximately constant and their values are given by 
([T8)) . We have also assumed that <SC tr so that the exponential can be approximated by unity in the 
integrand. From the above expression, we see that the necessary condition for the heavy modes to dominate 
at the initial time piitmi) < Ph{tini) is simply a > Hini (for 3 — 26 > 0) because this will guarantee that 
pi ^ Pini — Phitini)- Now, if wc usc that at high redshifts we have H ~ Hq^/Ti^{\ + Zmi)'^, with 57^ the 
density parameter of radiation, we find that the following bound must be fulfilled 

a » HoV^ril + ^«m)' ^ 3 X 10-36(1 + cV. (23) 



On the other hand, from ((22|) we easily obtain the effective equation of state parameter for the light 
modes. We see that they give a contribution to the energy density evolving as 

pi oc i^''-^ cx (24) 

where we have used that a oc fP. Then, the effective equation of state of the light modes contribution can 
be straightforwardly computed and we obtain 

wi=-^b+--l (25) 
3p p 
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Figure 2: In the left panel we show the cosmological evolution for the tower (solid blue line) where we can see that it tracks 
the matter behaviour in the early universe (we show radiation in short-dashed red and matter in long-dashed green ) and it 
exits that regime to eventually become the dominant component. In the right panel we plot the asymptotic equation of state 
of the tower as a function of 6 in a matter (solid line) and a radiation (dashed line) dominated universes respectively. We can 
identify the range of values of b for which we have the effective equation of state derived analytically in the main text in Eqs. 
12611 and | [27t . 

It is interesting to notice that, irrespectively of the background expansion determined by p, we obtain that 
for 6 — )■ 3/2 (which is our limiting case), the equation of state parameter goes to — > — 1. Of course, the 
result that we have obtained here is only valid for values of b such that wi < 0, i.e., for & > |(1 —p) (assuming 
positive p) because for smaller values of b the light modes will also contribute as a matter component and 
its energy density cannot decay slower than that. In particular, for radiation (p = 1/2) we obtain 6 > 3/4, 
whereas for matter {p = 2/3) the condition reads b > 1/2. 

In the right panel of Fig. [5] we show the equation of state for the light modes (or, equivalently, the 
asymptotic equation of state parameter of the full tower) obtained by using the analytical solutions given in 
(fTS)) and (fT6)) and, then, computing numerically the integral to obtain the corresponding energy density. In 
that figure we confirm the results found from our analytical estimates previously discussed. In particular, 
the linear dependence of wi on the spectral index given in (j25p that for a radiation dominated universe is 

wi\p^i/2 = -^b+l (26) 
whereas, in a matter dominated universe, it reads 

wi\p=2/3 = I - b. (27) 

We also confirm the breakdown of these equation of state when they become larger than 0, i.e., for b smaller 
than 1/2 and 3/4 in a matter and radiation dominated universes respectively. Thus, in order to have 
acceleration at late times wi < —2/3 we need b > 7/6. We shall see in the next Section, that this condition 
is not really necessary to have accelerated solutions since, once the tower dominates, the effective equation 
of state of the tower will asymptotically approach -1. 

From the previous discussion of the light and heavy modes evolution, we note that the total tower energy 
density can be approximately written as 




where we have used (|2ip and (j22p . Since we are assuming that Hint ^ cr, the dominant contribution at 
early times is given by the first term in the brackets, i.e., a matter-like behaviour. If we want the tower to 
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Figure 3: The 68.27%, 95% and 99% confidence regions obtained from the SNIa (blue), CMB distance priors (orange) from 
WMAP7 and BAO (blue). The combined fit is shown by the solid black contours. The fit has been performed by assuming 
that the light modes behave like a perfect fiuid with the equation of state given in I I27II and that dark matter is some additional 
field, not the heavy modes of the tower. In the right panel we show the vales of A and cr such that the light modes of the tower 
give the correct amount of dark energy today, given in Eq. II30I I. We have used fijv/ = 0-25 and b = 1.2, 1.3, 1.4 (solid, dashed 
and dotted respectively). The curves are not very sensitive to this value, being the major change due to b. 



dominate at late times, we need the second term to grow with respect to the first one. For a power law 
expansion a cxt^, this impHes b > |(1 — p) which gives b > 



^ for a radiation dominated universe and b > ^ 



for a matter dominated universe. On the other hand, by equahng both terms in the previous expression we 
can obtain the time t* at which the hght modes contribution becomes dominant, which is given by 



Hi, 



Hi, 



(29) 



This analytical expression has also been confirmed numerically from the exact analytical solutions. In 
the case of a cosmological constant, dark energy becomes the dominant component of the universe when 
H = ^/2^1aHo. This gives us the approximate bound H^, > ^/Wf^Ho for the light modes to dominate over 
the heavy ones early enough in the universe history. 

If the light modes are to play the role of dark energy, then we need its relative abundance to roughly 
accomplish the dark energy density parameter today, i.e, the following constraint must be satisfied 

1-26 



8ttG 

ml 



Pi 



3V2^(3 - 2b) 



Hj, 
a 



DE 



(30) 



This expression gives us the region in the parameter space in which the light modes of the tower can account 
for the current phase of accelerated expansion (see Fig. 

Since, the effective equation of state of the light modes only depends on the parameter b we can constrain 
its value from cosmological observations. We use the distance priors method as described in [17l | to confront 
to CMB data [Hi, the Union 2.1 SNIa compilation and the 6 BAO measurements given by 6dF [l^, SDSS 
21 1 and WiggleZ [22|- We show the corresponding confidence regions in Fig. [31 To obtain such contours we 
have fitted the mentioned cosmological observations for a dark energy model with constant equation of state 
depending on 6 as given in (j27p and assuming spatially flat sections. We show the confidence regions in the 
b — Oj\f plane, which then determine the remaining parameters A and a through the relation fij* = 1 — 
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4. Tower dominated universe 



In the previous Section we have shown how the tower can give a matter hke contribution from the heavy 
modes at early times and then the hght modes can take over and give a dark energy contribution. It is 
interesting to note that the tower could play the role of both dark matter and dark energy. If the energy 
density associated with the heavy modes is initially (in the radiation dominated epoch) dominant over that 
of the light modes, then those modes will behave like dark matter. Notice that under the aforementioned 
assumptions, a matter dominated universe is a consistent cosmological solution. This is so because, as we 
have seen in the previous Section, the heavy modes behave like matter in a matter dominated universe. This 
is not very surprising, since it is well-known that one massive scalar field oscillating around the minimum of 
the potential gives rise to a matter dominated universe. Thus, a collection of fields in such a regime, will also 
lead to a matter- like evolution. A difference between the two cases is that, in the case of one scalar field, the 
energy density oscillates and it is its average that redshifts as (1 -I- 2; )'^. However, in the case of the tower, 
the averaging is effectively performed by the integration over the different masses so that no oscillations in 
the energy density evolution appear. 

If we impose all the dark matter content of the universe be due to these heavy modes, then we need 

_ SttG i-wGa^A^ 

" 3i/o^V2^(3-26)(l + z„„)3 - (31) 

The tower will give a matter behaviour until the energy density of the modes that were initially frozen 
are relevant and, eventually, become the dominant component. Unlike in the previous Section, where the 
background evolution was assumed to be driven by some dominant barotropic fluid, in a tower dominated 
universe, when the light modes start dominating, its effective equation of state will tend towards -1. This 
can however been understood from our findings in the previous Section. While the heavy modes dominate, 
the light modes contribution will behave like a perfect fluid with equation of state parameter as given in 
(P7|) . When the light modes of the tower start dominating the energy content of the universe, the Hubble 
expansion rate decays slower (since the dominant component, the light modes, has a smaller equation of 
state parameter) and this will make the lighter mass components start decaying later. This will make the 
energy density of the tower to decay slower and slower, so that the Hubble expansion rate will decrease 
slower and slower. This process will go on until both processes are in equilibrium, which will happen when 
the de Sitter evolution is reached. 

Finding analytical solutions in the case when the Universe is dominated by the tower is much more 
involved than the case when the tower is subdominant. Thus, a numerical approach is required. To numeri- 
cally solve the equations, we discretise the tower and solve the coupled system of Einstein and Klein-Gordon 
equations for the associated discrete tower. Of course, the final solution will not depend on the discretisation 
procedure, provided the partition is fine enough. We have checked that this is indeed the case for our case 
and; moreover, we can recover our analytical findings in the previous Section when the tower is subdominant 
from our discretized associated problem. In Fig. 2] we show the obtained numerical solution for a universe 
dominated by the tower. We have chosen the initial conditions and model parameters so that the heavy 
modes are initially dominant and, thus, the tower behaves like a matter fluid. This is what we expected 
from our discussion above and has now been confirmed from our numerical solution. In this regime, the 
light modes behave like a perfect fluid with constant equation of state given by (P7)) . Then, for b > 1/2 we 
have that the effective w for the light modes is smaller than that of the heavy modes and, therefore, their 
contribution will become dominant. At this point, the mechanism explained in the previous paragraph starts 
working and the de Sitter universe is eventually reached. We can see this in Fig. |4l where the energy density 
scales like initially and then it decays slower and slower, eventually becoming constant. In reality, the 
asymptotic behaviour is a quasi de Sitter expansion rather than a pure de Sitter. The same effect can be 
seen from the equation of state that is initially and asymptotically approaches -1. 
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Figure 4: In the left panel we plot the tower energy density in a universe dominated by the tower. The parameters are such 
that the heavy modes dominate initially (6 = 1, crtini = 10), so we can see that the energy density decays like that of a matter 
fluid p oc (dashed line). At late times, the light modes take over and they asymptotically approach a (quasi) de Sitter 
regime. In the right panel we show the corresponding equation of state parameter where we see the initial matter-like regime 
(ui ~ 0) and the late time quasi de Sitter evolution with to ~ — 1. 

5. Discussion 

In this work we have studied a cosmological model with a continuous tower of scalar fields. Rather than 
obtaining the model from a particular theory leading to a continuous tower, we adopt a phenomenological 
model-independent approach and have simply assumed a Universe with a continuous tower of massive scalar 
fields, whose mass spectrum is given by a Gaussian distribution peaking at the origin and with a certain 
dispersion a. For the cosmological evolution of the tower, we have used initial conditions according to a 
power law with amplitude A and spectral index b. Then, the behaviour of the tower during the different 
epochs of the universe expansion history has been studied. We have found that the heavier modes of the 
tower give a matter-like contribution, as expected from the usual results of quintessence modes. However, 
for the light modes, we do not obtain the same behaviour as for one single light scalar field, that would give 
a A-like contribution today, provided its mass is < Hq. Instead, we obtain that those modes cooperate to 
behave as a perfect fluid whose equation of state depends on the initial spectral index b and ranges from 
to —1. Only when t» — |, the light modes behave nearly like a cosmological constant. We thus avoid having 
to introduce a very small mass like in the usual quintessence models in order to have acceleration today. 
The reason for this is that we have mass modes infinitely light in the tower that will remain frozen and, 
then, will give a constant energy density contribution, no matter how small the Hubble expansion rate is. 
In fact, this is the mechanism behind the effective equation of state parameter. 

We have obtained constraints on the parameters of the model from CMB, BAO and SNIa observations 
by modelling the light modes with a perfect fluid with constant equation of state. The tightest constraint is 
obtained for the spectral index 6, while the initial amplitude and the distribution dispersion are degenerate. 
The reason is that these two quantities only contribute to the density parameter today of the light modes, 
while the effective equation of state is entirely determined by b. 

Finally, we have studied the case of a universe dominated by the tower. In the case when the heavy 
modes dominate, the universe behaves as dominated by a matter-like component, as expected. However, 
for appropriate values of the spectral index b, the light modes will eventually dominate and their effective 
equation of state approaches -1 so that the universe eventually becomes de Sitter. 

We can also comment on the fact that having a continuous tower with all possible masses, including a 
massless component might be troublesome. If a massless (or very light) scalar field is coupled to standard 
model particles, then it will mediate a long range force that has never been detected in fifth force experi- 
ments. Ways out to this are usually given by making the field dynamics environment dependent so that its 
effective potential (or mass) screens the field effects in appropriate environments (chameleon or symmetron 
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mechanism) or non- linear terms beeome important in high density regions (Vainhstein mechanism). These 
mechanisms could be implemented for the individual modes of our tower to avoid fifth force constraints. 
Another potential way of avoiding such tests could be developed by taking advantage of the distribution 
weighting the different mass modes. Thus, one could try to make the mass distribution environment depen- 
dent such that it strongly suppresses the lightest modes in regions of high density and, then, no fifth forces 
will arise in local gravity tests from the tower. 
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